5 


i)Clc 


THERMOMECHANICAL STUDY OF COLD ROLLING 




B . R . VI NAY KUMAR 


Ih . 

KUM 


L. 



department of mechanical engineering 

INDIAN INSTITUTE OF TECHNOUOGY, KANPUR 

JUNE, 1988 



THERMOMECHANICAL STUDY OF COLD ROLLING 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

Master of technology 


feroeox 


by 

B . R . VI NAY KUMAR 


CO the 

DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

JUNE, 1988 



jViE 



JAN 1989 


No. 


103071 





jt 


l^iSiO- (V\ ^the 



^ !L» 



This is to certify that the present work on 
"Thernomechanical Study of Coir' I’lOlling ” by S.R.Vinay i^uiaar 
has been carried out under our suoervision and has not been 
submitted elsewii ^re for the av'ard of a degree. 



(A'l.K.i.luju) 
Professor 
Mechanical i£ngg. 
I.I.T. Kanour. 



- V^e rv. -- -(‘V" 

( P.M. Dixit) 

Assistant Professor 
Mecnenical Engg. Deot. 
I.I.T. Kanpur. 



ACKM0HLE33Evl£N TS 


I express my sincere gratitude to Dr.M.K. Muju 
and Dr. P.M. Dixit for their constant encouragement, 
invaluable guidance, comments and criticism throughout 
the present work. 

I am thankful to Dr. S.G. ISiande for allowing the 
use of CAD-project facilities. Thanks are also due to 
Dr. T. Sundararajan and Dr. K.K. Saxena for their valuable 
suggestions « 

I also thank Ramdass, CCS, KGS, Mrs. Biswas, Prakash, 
Satish, Nigam and Lekhraj for helping me tide over technical 
and non-technical problems in CAD-P. Thanks to Mr. Sen, 

System Supervisor CAD-P, for the delay in the thesis. Thanks 
also to ND-560 for standing through thick and thin with me. 

I also wish to thank my friends Diva, Venki, Koppaka, Per 
Sambu, Raju, Karra, Venu, Pokuri and Rampa for making my 
stay at IIT - Kanpur a pleasant and memorable experience. 

Finally thanks are due to Mr. U.S. Mishra for his 
neat and accurate typing. 


f 


I 


June, 1988* 


-B.R. VINAY KUMAR 



co'fThiir.o 


Paqs 


TJGT 01- ^lOORHS 
MOf'AF.MCLATf RH 
ABSTRACT 


CHAPTER I IHTPnDiTo-jrTOr-I 

Prcivlous W'^rk 
Scope of Pronrint V/ork 

C!-IAPTER II HA THEM ATI CAL i’DUE LLIHO OP ROLLIHCi PROCESS 

2»t Mechanics of Rol.llnq Procoso 

2.2 Covorninq Equations 

2.3 E ("lu a t i on s for S toad y S c t e T \-io 
Diniensi on.al PcobJ.en; 

2.4 Penalty Formula lion 

2.5 Poundary GoncJltiooc 

2.5.1 Entr-y and Exit noimdarles 

2.5.2 Convective ''boundary Sides 

2.5.3 i'.ol 1 “Contact boundary 

2.5. 4 Ax is of Sy mifi '■ try 

2.6 Hondiirtons loncoLioat i on 

2.7 Galeri-in ’^^ornutlalion 

chapter TIJ FTHTTE ELFMEr-’T MFDrHXJNn APD IMPLEMENTATION 

3.1 D;i screi.'i notion ami Shape Functions 

3.2 Elmnontnl Area F.xproscionc 

3.3 01o!.‘'al Assembly of Fl^menta.l Area 
Expression 

3.4 Eloinent Boundary Expression 

3.5 Evaluation of and _fy on the 

Boundaries ’™“ 

3*6 Global Assemljly of Rlcrnental Boundary 
Express 1 ons 

3.7 Global Equation 

3.8 Boundary Con'iitions 

3.8.1 Special Boundary Con'lihions 

3.8.2 Essential Boundary Conditions 

3 .9 Programme Implemontation 


1 

2 

4 

5 

5 

7 

1 1 

12 

13 

13 

15 

16 

17 

18 
20 

25 

25 

29 

35 


42 

45 

46 
46 

46 

47 
47 


CHAPTER IV TEMPER ATiiRE MFA3HREMFNT TM ROLLING 


P.O 

Or, 


CHAPrm V RESULTS AMD DIGCiJSSJUid 


REFERENCES 



LIST 0^ FIGURES 

Fi.gure Title Page 

2*1 Distribution of velocities in rolling 23 

2.2 Rectargular Cartesian coordinates 'for the 24 

domain 

2.3 Boundary codes for the domain 24 

3.1 Typical element e and the local node numbers 49 

3.2 Typical boundary element b and the bbundary 49' 

local node numbers 

3.3 Backv/ard and forward slip regions on the 50 

roll-v;ork interface 

4.1 Arrangement of thermocouples on the workpiece 52 

4.2 Schematic diagram of experimental setup 52 

5.1 The 56-element finite element mesh ' 60 

5.2a Variation of V along the Interface 61 

Reduction ratio = 15.38f/^ 

Angular velocity of the roll = 2.0944 Rad/sec. 

5.2b Variation of temperature along the interface 62 

Reduction Ratio = 15.385^ 

Angular velocity of the roll =2.0944 Rad/sec. 

5.2c Experimentally recorded graph reduction = 15.38?b 63 

5o3a Variation of Vg along the interface 64 

Reduction ratio = 27 % 

Angular velocity of the roll = 2.4 Rad/sec. 

5.3b Variation in Temperature along the Interface 65 

Reduction Ratio = 25/o 

Angular velocity of the roll =2.4 Rad/sec. 

5«3c Experimentally recorded graph Reduction = 25.00% 65a 

5o4a Variation of V along the interface 66 

Reduction ratio = 21 .86?4 

Angular velocity of the roll = 1.92308 Rad/sec. 

5.4b Variation of temperature along the interface 67 

Reduction ratio = 21 .86% 

Angular velocity of the roll = 1.92308 Rad/sec. 



Figure 


Title 


Page 


t.5a Variation of V along the interface 68 

Reduction ratio = 21 .86°o 
Angular velocity of the roll = 2 <.40384 
Rad/ sec . 

5.5b Variation of temperature along the interface 69 

Reduction ratio =21 .86^ 

Angular velocity of the roll = 2.40384 rad/sec. 

5«6a Variation of V along the interface 70 

Reduction ratio = 29 *40^3 

Angular velocity of the roll = 3.8461 Rad/sec. 

5.6b Variation of temperature along the interface 71 

Reduction ratio = 21 .86/o 

Angular velocity of the roll = 3.8461 Rad/sec. 

5o7a Variation of V along the interface 72 

Reduction ratio = 29.40" 

Angular velocity of the roll = 1.92308 rad/sec. 

5.7b Variation of temperature along the interface 73 

R-eduction ratio = 29 .40?o 

Angular velocity of the roll = 1.92308 Rad/sec. 

5.8a Variation of V along the interface 74 

Reduction ratio = 29.409^ 

Angular velocity of the roll = 2.40384 Rad/sec. 

5.8b Variation of temperature along the interface 75 

Reduction ratio = 29.40“o 

Anaular velocity of the roll = 2.40384 Rad/sec. 

5o9a Variation of V along the interface 76 

Reduction ratio = 29.4093 

Angular velocity of the roll = 3.8461 Rad/sec. 

5.9b Variation of temperature along the interface 77 

Reduction ratio = 29.40^3 

Angular velocity of the roll = 3.8461 Rad/sec. 

5.10a Variation of V along the interface 78 

Reduction ratio = 33.75% 

Angular velocity of the roll = 1 .92308 Rad/sec. 

5.10b Variation of temperature along the interface 79 

Reduction ratio = 33.75?-o 

Angular velocity of the roll = 1 .923 0 8 Rad/sec. 



Figure Title Page 


5.11a Variation of V along the interface 80 

Reduction ratio = 33.75“i 

Angular velocity of the roll = 2.40384 Rad/sec. 

5.1 1b Variation of temperature along the interface 81 

Reduction ratio = 33 .75Pj 

Angular velocity of the roll = 2.40384 Rad/sec. 

5o12a Variation of V along the interface 82 

Reduction ratio = 33 .75°^ 

Angular velocity of the roll = 3.8461 Rad/sec. 

5.12b Variation of temperature along the interface 83 

Reduction ratio ^ 33.75fo 

Angular velocity of the roll = 3.8461 Rad'^sec. 

5.13a Variation of Vg along the interface 84 

Reduction ratio = 16'o 

Angular velocity of the roll = 2.015748 Rad/sec. 

5.13b Variation of temperature along the interface 85 

Reduction ratio = 16°o 

Angular velocity of the roll = 2. 01 5847 Rad/sec. 

5.14a Variation of V along the interface 86 

Reduction ratio = 

Angular velocity of the roll = 4.013496 Rad/sec. 

5.14b Variation of temperature along the interface 87 

Reduction ratio = 16^^ 

Angular velocity of the roll = 4.013495 Rad/sec. 

5.15a Variation of V along the interface 88 

Reduction ratio = 16Pi 

Angular velocity of the roll = 5.03937 Rad/sec. 

5.15b Variation of temperature along the interface 89 

Reduction ratio = 16°o 

Angular velocity of the roll = 5.03937 Rad/sec. 

5.16a Variation of V along the interface 90 

Reduction ratio = 24 % 

Angular velocity of the roll = 2.015748 Rad/sec. 

5.16b Variation of temperature along the interface 91 

Reduction ratio = 24 % 

Angular velocity of the roll = 2.015748 Rad/sec. 



Ficpjre 

5.17a 

5.17b 


Title Page 


Variation of Vg along the interface 92 

Reduction ratio = 24"'^ 

Angular velocity of the roll = 4.031496 
Rad/ sec . 

Variation of temperature along the 93 

interf ace 

Reduction ratio = 24^ 

Angular velocity of the roll = 4.031496 
Rad/sec. 



MOMBICLATIIRE 


A 

4 

B 

B 


c 

e 


hl 

K 

lb 

N 

~P 

LIV 


£ 

q 

R 



Assembled global matrix 

Unknown pa?'t oF RHS containing and tI* 

Assembled riqi>t hand side 
Mon-dimensiona Used strain matrix 
Specific heat (J/kg- C) ’ 

Internal energy (j/kg) 

• Initial half thickness of the sheet metal (ra) 
Final half thickness of the sheet metal (m) 
Elemental area matrix 
Thermal conductivity (¥/m~°C) 

Length of bha boundary element 
Shape function vector for temperature 
Shape function vector for pressure 
Shape function matrix for velocities 
Hydrostatic component of stress tensor 
Elemental nodal pressure vector 
Heat flux (W/ra^) 

2 

Normal heat f^ux (W/m ) 

Reduction ratio 
Radixis of the roll (m-) 

Temperature ( 

Elemental nodal temperature vector 
Normal traction on the boundary (N/m^) 
Tangential traction on the boundary (N/ra^.) . 
x-direction traction on the boundary (N/m^) 
y-diXGction traction on the boundary (N/m^) 

\ 



TRACT 

'•"n tfip pt'o'^onl: Iho'^in pn p! f mode to pr<?d 5 ,ct. 

the pressTjre, velocity -'Mfj "cc:pe:cptuTe di c- hribution in a 
ntrip dorirsa cold r'^l.'’<i-nq. Tfio c'^nt Inuityj v.omentum and 
energy t’qtipt .ion'^ wore u'^od Oor r c.rrriiX riino the problem 
theorotj CP 1 ly« Th'^ Frict'on on tJip roll-'./ork, interface wac 
assumed to be oF siicki.nq ru-;n’ slippinq type* The 

mechanic a.l and thermal p->'rtp''v i-,-? na or v/orlc v/ere assnmou to 

be temp'^rp tm:e dop'-n 'I'-nt • Tine h'''p!', was assumed to dissipate I'V 
convect''on on free socles oikI by r.u-i.juction on roll-work 
interface. The sot of non-lln'^ar equations resulting From 
the above formulrati on were then soJvofi iteratively using 
FEM proco' lure based en Gplerhin's nporoach until conveP'genco . 
The validity of the thooretical medel is coifirmed by comparing? 
the temperatures exp<^r1inen tali 7 . The velocity profiles also 
turn out on expected linaso 
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aiAPTHR I 


IMTnOD‘T''Rr.TON 


Mumerical and computer aided analysis of forming 
processes have recently started to make an important 
contribution in industries. From the physical point of 
viou, the analysis of forming processes poses a range 
o" complex problems in continuum mechanics and heat 
transfer. Solution of the problems involve usage of 
sophisticated numerical techniques. The finite element 
method eme.rges as the preferred approach and its flexibility 
leads to the development of powerful programme packages 
capable of treating wd de range of problems. 

The forming of many indtJistrial products can be 
accomplished by many different process all of which 
involve large, permanent deformation. Rolling is one of 
the most important industrial process in this respect. 
Rolling may be done as hot or cold rolling. Rolling 
operations done above the recrystallization temperature 
of the material are said to be hot rolling process and cold 
rolling process if done below recrystallization temperature. 

So far the traditional approach to predict the flov/, 
roll pressures, forces and torques is to establish an 
average temperature in the roll gap and use it in the 
determination of an nvorngo flow strength. Since this 



mothod of nnn.tysi? lines not tote into nccount the heat 
generated .dn<^ to elthnc filast5.c de-'^oTipati on or friction 
at the intorfacG, the depondonco of yiold strength and 
other thermally senaltlvo proper i; j eo of the material 
cannot be ooti of oc to’’"iiy incorporntori, in hhe analysis* 

The thermal gradients in theinatoriol are dependent on 
the roll-material used, reduction ratio, angular velocity 
of rolls etc* 

Previous Hork 

The importance of tomporaturo calculation have 
been long recognised., ^titil recently, the majority of 
the vrork has been to dvocouple the problem of heat transfer 
from the metal defoT*mation nroblera* The approach has i>onn 
to determine the flov/ fields in the problem either crpe.ri- 
mentally or by calculations and then by using these fields 
for the calculation of heat goneratloru 

From earlier v.^orks [1] on metal~-So rndng analysis, 
two different approaches have emerged: the ’flovd approach 
and the * solid approach In the flov; approach, the metal 
is considorod to behave as a non-Idowfconian fluid and the 
finite element eguations are solved accordingly. This 
m.ethod, no doubt being very effective does not take many 
subtleties of the el astic-plnstie constitutive equations 
into account and is not in a position to predict residual 



stresses. Tn the solid approach, the materiel is considerrri 
to behave as a classical elastic-plastic solid and the 
finite element eqnation^- are solv'^d v;ith methods used in 
small-strain plasticity prnhloms® 

In 19785 Allan, 8hah and Lnholi did a computes.- 
aided simulation [ 2 ] of strip rollj nq process by combining 
analysis o" rolling m'^chnnl cs and by developing a finite- 
difference method for predicting time-deponden t temperatures 
in the strip and the work rolls. In 1981, Zien Kiev/iez, 

Onate and Hein.rich car.Tiod out a coupled thermomechanical 
analysis [3] for flow of metals using flov/ formulation® Pr'^sma 
thesis has used the idea of flow formulation from this nap--^c,.'q 
notable v/ork in this area has been carried out by Cornfiold 
and Johnson [ d] , in hot rollinq, Johnson and Kudo [ 5 ] using 
upper-bound solutioji techni que® On the (^experimental ■ side. 
Bishop, Piice, Jesvn.et, Altan, Karagio/J.s, I.encard and others 
have carried out various experiments in determining the 
temperature distribution and finding out .roll-interface 
temperatu.rGS . Jesev;iet a'ld Rice [6] had embedded thermo- 
couples into the v/oric rolls and evaluated the temperature 
distribution on the interface* Knraqiozu cind Lenard [7 ] also 
used em.bedded thermocouple techniquo but they emibeddod 
it in the x^orkpiece to evaluate the temperature di strii->u tions * 



Scope of F'-rcoo nt \!oT k ’ 

In the preset! h '•'^rlc, the nrii.n thi.mst has been to' 
GvDluote the velocity or flovr uistr'i butj.on and the 
temperattire d'l stribntion in cold rollinq process® Pre’'?-!- 
ously mO'':t the r ''’-a o archers who have worked in the nt’mer-- 
ical thorrnoincehani cal onnlys.-r s of cold rolling had used 
the techniqno O'^ docampling the itiochanics and heat transfer 
problot’tCo /hi onk I ov;] cz e al®, fm' tiie first time appl;* c'’' 
the technique of cottpling thornial and mechanical problo'ijc 
ttsing flow rormnla/tlon » or the .s'"-'<''^5T'o 1 arsuraptions they 
assumed a *stjc!:lng' condition on tfv’ roll- Interf ace® 

This implied that the velocity of t'le v;ork material on 
the interface was same as the ro.i l-sucface speed. The 
present thesis differs from Zienkl owi cz hs v/ork in that 
slipping is assumed to occur at the interface and Coulorr.bic 
frictional forces v;ith some upper liiriit as the sticking 
friction condition® Further more, the frictional constant 
is asstmed to be constant over the interFace due to lac!: 
of exact data* 

The neutral zone v/as assumed to be a point® The 
mode of heat transfer on the free surface is assumf'd to 
be through free or natural convection® Ti)e experiments 
which were carried out in the present vjork viere aimed rot 
measuring the interface temperaturoo 



CHAPTUR TI 


fMT HR'.'iATTC AL M*^D'-LLI * p j^^yrT. .; ppoCH'^S 

Tn ’.hip chapfcor tho ic.odol "ot tho rolling procero io 
developod. Tht^ nrocosc ’ o rnodoilori op a .Pteady-.state io-'o 
dimon.si nnal na^oblern® Tb. o interaction vjith the roll and 
the .en.rr- 'ondi ng.s in x-oor e.ocnted by appropriate boundary 
conditions. The con.sti tutive onnahlen a Tor mechanical 
and th<^rmal hohavioux- of the ohrlp n.y'a si. abed. Mon“ 
dime n-cionpJ i nation oP tho govorni no oguation.a and the i>?oai: 
formulation is pre.sontod at the on'i oP this chapter. 

2-1 Mo c ! t a nips of R o .7. l.ln g Proc o sjs 

Many experimental and theoretical invest igation.s 
have shovm that the volocltieo, oti'nin-s, stresse.s assoc i a t-^d 
with motion of a rolled ,strip are not uniPormly distributed 
across the sections. 

The material being rolled unriergoes varied velocity 
and strain distribution (see fig. 2«.1)« The velocity across 
the section of the strip in the unstrained zone is distr? - 
'buted uniPormly, whilst in the contact zone of deformation, 
the outer layers of the metal touching the rolls terxi to 
move faster than th^ inner lay-^r.So The outer layers of 
the metal in contact veth the roll have less velocity than 



the roll nf^ar the enti'y. Put the velocity of tho outer 
layers (denotf^d by v^) cecreases a e it reaches the exit 
and finally h'-corses more than the spef^d of the roll 
(denoted by v^,) nea'^ tto exit- Tho zone of transition 

uhere v becomes equal to v^ is called tho nei-itral zone* 

s r 

Even thouqh experimental evidence shovrs that this transition 
occurs over a zone, many researchers have considered tills 
zone as a point in their analysis Tot' Jacl: of exact data 

and the difficulty .in meriollinq i i, as a. zone* The zone 

in v/hich v„ is les.s than v is call'-rj bnokvrard slip renlon 

and the Prictional force on tiie iiiotaJ is in the direction 

of rolling. The zono in vdiicl) v,, is aroater than v^ is 
called tho zone of forward slip roqion and the frictional 
force acting on tho strip t-ri os to retard the motion. 

Many researchers have come witri an analytical exprossi 
for determining the neutral point. Notable among thorn and 
given below are? 


(1) Avitzur’s formula [8] t 
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( 2 ) Chaturvneii ’ s Pci'mula [9. 
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v/here, 

'a " '’o 

r = rodHct5.on rat/i^'i -■■ — *r" — '■™"* 

■ V> 

= .Initial tlnichnnaa ata'i p 
= final t'licknosa of atrip 
= non+.ral amlo 

‘fi 

'? = friction p.nale 


& = angle of contact 

R = radio a of the roll 


m = f.ri.ction factor 



back tension 


= front pull® 


Avitzur’s formula v/ac chosen to calculate entry and 
exit velocities of tho strip. The choice was based because 
this particular formula was most qonoral in form and 


2 .2 Govorni.nq Equations 

Tho thorir.ornocha n1 c >1 hofiovinin' of a material is 
governed by the following throe onontionsr 



(1) Connorvation of mnor; ov continuity equations 

P 'I P div'' 7 ~ 0, (2 •4) 

( 2 ) The bnlanco or rnr'i.i'-'nturn or oquotion of motions 

♦ 

Py — V e ^ ~ ^ ( rov zero body forces) C2®f 

('^' ) i4ons'"0'v-'tion of enori'iy c' thr- •rifct law of 

thori-odynnmic c ; 

P o -! V • q •" .? » V “ 0 ( « f'' ) 

Bosidos the abovo equations, v/o ulao need the cons t i.t ui iv o 
emulations « The strip briioviour is modelled by the folloviina 
constitutive equations t 

( 1 ) Stress-strain rate relation! 

In a flo\^ process, the def^ormation is normally repre- 
sented by the strain-rate tensor (ncQlocting elastic 
strain) 


e = ^ ( vY_ + V v"^ ) 

. ’ 

Let B rond e be the deviatoric po^ts of the stress and 
strain rate tensor! 

^ ' 1 

a = -Pi-! 5 , p “ tr 

^ 4 ^ 

g = (“ tr |> 1 -!- £ 


f 



Then they ore hy 


• ! 




where the vir,cos5.ty f-.l:' given by 

1 

}i ^ ~7 

"or fi rlglo' f>lor,tle n’. al v/hJ ch /iel'.in according to 

Von-iMlGCC c’^lt or I on. 

( 2 ) The dennl ty"-tf'ir.per<?.t'.r''e ro.l.'itU>ns 

In a plactic inatorlal, vo.lomo ahnncfos only v/hon 
there is n change in tempornturo To Th^rf^-fore depf’nde 
only on T« 

T 

P = P exp [ - / 3a(s) sic ] (2.2) 

'^o 

where « = cne<^ficient of linear expansion of the metal 

Pq = density at room temperature 

T„ = room ternoerature. 

o 


( tr ^ ) 


2 



(3) Meat •Plii>:--tempo.Tnt”'’‘e aradl'^n !• falatien* 

2 = -~r:vT (2.9) 

wh ere , 

K - thermal condnctivlty- 


( 4 ) Specific internal energy-temper a ture relation: 


T 

e = / c(s) ds ■ (2 -1 O' 

0 


v/here» 


c ~ specific heat* 

H-^-re the material properMae a , p. y a and c are 

y 

cnnsiderod to bo tomperaturo d'^p^nd.enfc ® 

tubstitnting the conetl tntive equations [ 2 . 7 - 2 .1 0 ] into 
the governing equations [ 2 . 4 - 2 . 6 ] ue got the follov;ing 
equations: 

+ lyr) r.- 0 (2.11) 

- V. 5-0 (2.«12> 


/3T 

V. V - 3 

P (|| + V . vV 


and 
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2*3 F.nu.Ttlon s ~^or G ' /> C'-Iy Ft ate Two Dirfion^ional ProblciTt 

Th ncy.'ip i a^si’rti^d t.o i'C vc^y broad in relation ’/o 
the lenqth of tho ny- of con inch and the thickness of the 
strip. ConsSariup nl- ly tho offoct or '^prr-ad may be neqlocted 
and tho nroblom can bo consi^-’oror! pc two-dimensional flat 
rolling. 


Roctnngnlar can’tosinn com'c’inato system as depict'^'"' 
in fig. 2.2 is conaiciorod, vri th the X~ar'S along the lonolh 
of the strip being rolled and I'-axis along the thickness of 
the strip. 

For a 2~D prohnloiri bf'io s/’olocity vector has only i;v '0 
components and honce 

V = u i -1- V j 


while tho stress and strain-rate tonsor have the folloidno 
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Tho governina ^atiationo f2»11 -- 2e13l are then 
reduced to; 
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2.4 Penalty Fornrmlation 

On observing equation [ 2.14] it is observed that it 
devoid of pressure terms. Because of this the finite element 
formulation of this equation results in a matrix with zero 
diagonal terms. If the number of zeroes in the diagonal are 
significantly large, tho coefficient matrix becomes singular# 
To overcome this difficulty, we modify tho equation [ 2«14 ] 


as follov7s; 
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where X ie a larne nn>'’onr. In the limit X oo 
entiation [2.17] reducaa to oq. [2.1'!] . 


2*5 Bonndary Conciitions 

Tho boimdnry conditions '^peci'ied on any part of 
the boundary can b<^ claGsirieci ns 

(l) osGonhial boundary conditions i.e., specification 
of teifiporatux'e T <in'’l ^/oloc' hy ''■£« 

( 2 ) natural boundary concii ti one i.e., spocifi cations 
of stress vector t - ? « n and normal heat 

flux - £[ * • 

( 3 ) mixed type nro compon^nt'^ O''' V and the remaininq of 
_t aT-e specified. Sometiin'^s thc' ^r'^lation bet\-/oon 
the components of V and _fc are np^cified. 

2.5.1 Entry and Exit Boundaries 

This is the roqion sh.ovm by sides 1 and 5 in the 
fig. 2.3. ThevSe boundaries are assumed to be far away 
from the zone of deformation. Thus the temperature on 
these sides can be assumed to be at rnoir. temperature. 
Furthermore the y-velocity component of velocity 
vector V) is zero on both the boundarj es. The X-voloclty 
components (u]' are calculated from the roll surface speed 



and AvitziJrr; formula [2»1 ] 


The n-volocity components are found as follows? 

First the neutral angle ® is calculated. Then the 
following continuity equation is used in obtaining 
u-velocity component at exit. 

V Co s 

^ [h^ -( 2F(1 - Cosa^) ] (2„10 

o 


Then the velocity at the entry is: 


u 


1 



In actual practice it would I.-^e more appropriate if 
the u-rvelocity components at both entry and exit are found 
experimentally. 

Thus to summarise, the boundary conditions are: 


T = T on both side 1 and side 5 
oo 


u 

V 


= u 


0 


on siae 5 


and 


u = u. 


on side 1 


V 


0 



I (Ion 


IT) 


^ *2 r"o n v r' c 1: 1 \r o n i *11 .» ( 1 . 1 r y 


The Gides nri v'h?ch conv^^c tvive boundary conditions are 
to be applied are the -'idoG 2 and 4 as shovm in fig* 2«>3» 

Tiae losses duo to radiat'ion are considered negligible® 

The heat that Is being convoebed away to atmosphere 
by free cnni.^-nction Is given iiy 


Cl = h ( ) 


( 2 .. 


where h “ heat transfer cocFf I cisenb given by the formula f 1h] 


/Alt 

h - 1 ,?2 (“T*) 

W 


0®25 


( 2 .: 


where , 

AT = temperature difference betv/een the free surface 
and room temperature, 
w = vridth of tlie plate. 

Since the strip velocity is along the free surface the 
following conditions are considered! 

V = 0 and ~ 0 

Then the boundary condition's can bo summarised as follows s 


V V, 0 

- 0 

q -h(T-T^) 


and 
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2.5.3 Roll-Contact Rounnarv 


This interface io shown by the side 3. On this 
surface, due to the frictional work, heat is generated. A 
part of this heat is conducted into the rolls. The rest 
goes into the strip. The expression for the heat conducted 
into the rolls is: 


q 


r 


- TJ 
211 


v/here, 

Kj, = roll conductivity and 
R = roll radius . 

Assuming that the diametrically opposite end of the roll 
to the interface is at room temperature. 

The heat generated by friction is: 


q = tgVg ; 

where , 

tg = frictional force and 

Vg = tangential velocity of strip. 

» 

Therefore the heat flowing into the metal is s 
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Hie mechanical boundary conditions for velocity and 
tractions are bit complex in the roll contact zone. The 
normal velocity component along the radius v^ is zero 

V = 0 « 

n 

Mixed type of friction model is considered ;daich 
incorporates both slipping and sticlcing friction. The 

relation between tangential force t„ and the normal force t 

s n 

i s given as 

t- = P-t„ if t < (slipping friction) ■ (2»7;'’s) 

s n s 3 

and 

t_ = — ^ if t > (sticking friction) (2<-2P,b) 

® 3 ® “■ 3 

is considered to bo constant along the arc length. 

In the sticking friction zone the layer of the metal in 

contact v/ith the metal moves vrlth same speed as the roll 

and the frictional foreo acting is equal to the yield 

shear stress. This type of formulations leads to unknoi^m 

terms in t^^ and t^ in the expression for on the 

interface. The elim:] naf I or\ of those unknown terms is 

dealt with in detail in the next chapter. : 

i 

2.5.4 Axis of Symmetry : 

The y-velocity component v is zero along the axis 
of symmetry, i«e., side 6. Furthermore, t^ is zero along i 

the side. The heat transfer across this boundary is also I 


on this side 


can bo ^ . 

-iii.n.arxsed as? 


V ~ 


"C, 


end 


n 


0 

0 

0 


2 . 6 Mon- cb.rr,on -i n n t 



f-^^n-dimonsionain. .. 

is usually rjone to avoid ill 

conditioning and nu,,.,.-, ,,, ,. 

,, , ' O ' I. f I cul hies . Tn this section n 

the physical quantity n. an 

--0 are non-dirr,onsi onalised using some 

charactori sf i r- a • 

' • • c ^ Mr< o n “I o n 1 

•nnj. q, related to the prooes 

parameter .Ijic, ^ 

roll velocity V , 

Sheet h , amb'lonf ^ 

o ni I.oir.pornture t 

DO 
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2 L 
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f:5nnl thickness of the 
o\:c<^ Tints 




V 


V/ 
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* T — -i* 

y i rp 
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The mechanical and the ' 

nondimensionallsed 

with rospoct to thoir ^ 

room teujpornture values i 


lr~* * iT r-: __li 
o p - 
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; R 
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'o 


K. 

k. 


C 2 »23 


^'fhere 



and , P » a » c and K are hhe room t omperature values* 

yo o o o o 


The presstjre and th^ stros*^ tensor are nondimensionalised 


by the relation; 


/®«*0 w.l? 


ix 

I o 






Substitutin'! oqusatlons [p«2Sa,2.2'’l'','--»^l’'Sl into equations 
[2.15a, 2«15b, 2.16, 2.17] v/o gel:? 


(e + e } _ ga X f cr (i 
^ xx vv^ o no i- ' 


] H- 




p- (5 ^ , V ^i-)] - -I 

^ o 3 X ay ^ ^ 


-XX.) 


5 a) 


rp (H ! V 

P Q 0>: 3 y 


^ ***' 3 CT 

3 ji)]_ ( Si 


25.,.y.y) 


(2 *255 )i 


P c T„ h 


o 0 00 0 


M- n r 


p C (u '! 


- 3T^ 

V — } 


K T 

0, . 00 

0 


[ --r (K 

3 X 


3 X\ — tt 
— i) -I p ( e. 


[ ^ XX 


X ( E 
3 XX 


(2j.2o ) 



The nondirnonnxonal normol heat flux and the stress 
vector are: 
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n 
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t rr 


t- 
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t 
w r 
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o 


V , 


2»'7 G a 1 r k -in F o rmu .1 a b x o n 


Thr^ '^ovorninq <^nunt1rins i.n Ih? iiiocJcl are rron-linonr 
an d varl at. i onnl Formn.latlon dooa nor oxint hare. Galerhln'a 
weighted residual mehhod is adopted in the present case. The 
weak form of the prohlorn is dori.ved in this section. The 
Finite Element Formulat.ion ot this folloxrs in the next 
chapter. 


Let u, V, p and T be the functions which satisfy 
the equations f2«2d - 2.26] and aJ.il the boundary conditions 
exactly. Then they v/ill con.stitutn a v.'onk .solution of 
the problem if the follov/ing integrals are zero! 
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vrhere Wp, and V!j are the vroiqht functions which 

satisfy the homoqfoneous versions of tho boundary conditions® 

Inteqrating by parts certain terms and using the 
relations 

_ I 

= -p 1 -1- 2ilp Q e’ ? 1 = ‘^•n and = Kv T . n 

and the boundary conditions on the weight, we get the 
following vreak form of the problem, 

dA = 0 {2®28) 

-'■p I 2 dA= ^ T 4 dA 1- / dA 


( a 29 ) 


/ lo dA 
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dA 




'C 


v;here 
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(2.31 ) 
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1 - the velocity of the outer portions of the strip cross section ; 

2 “the velocity of the middle portons of the strip cross section ' 

3 “the mean velocity of the strip cross section; 4-veloclt dia“ 
gram for the unstrained zone • 5 -velocity diagram for the defor- 
mation zone at the entery, av^ay from the contact zone; 6- velocity 
diagram for the zone of backward slip; 7- velocity diagram for 
the natural zone ; 8- velocity diagram for the zone of forward slip; 
9- velocity diagram for the deformatron zone at the exit, away 
from the contact zone; 10-velocity diagram for the unstrained 
zone at the exit. 


Fig. 2.1 Distribution of velocities in rolling. 



R 



0 

Fig. 2.2 Rectangular Cartesion coordinates for the 


domain 



Side 6 


Side 1 - Entry Side 4 -Convective boundary 

Side 2 - Convective boundary Side 5 - Exit 

Side 3 - Roll - work interface Side 6 - Axis of symmetry . 

Fig. 2.3 Boundary codes for the domain. 


CHAPTER ITT 


FINITE ELEMEMT MODELLTMC AND IMPLEMl-MTATTON 


In this chaptox', the method of discretisation, shape 
functions and the elemental area expressions are discussed. 
The method adopted in applying the essential and the 
natural boundary conditions especin.lly the conditions on 
the roll-vrork interface are discussed. 

3 .1 Discretisation and Shape Functions 

We use the convergence of the variational FEM in 
selecting the shape functions or the conv^ergence criteria 
for the Galerkin FEM are not well established. The order 
of the highest derivative of p, {j, v and T is respect- 
ively zero, one, and one. For the coir.pletoness criteria, 

Pf V, T and the first derivative of u, v and T should 
tend to pon-zero values as the element size goes to zero. 

The compatability criteria places no restriction on p 
but demands that u, v and T be continuous across the 
interelement boundaries. To satisfy the above conditions, 
a zeroth degree polynomial for p and linear polynomial 
for u# V and T suffice. However, zeroth (iegree polynomiol 
is not feasible in FEM. Therefore we use a linear polynomial 
for p. To avoid certain numerical difficulties, the 
polynomials for il, v and T should he one degree higher 
than that for p. Therefore we choose quadratic polynomials 



■For ii, V and T. For thio v/e nood d-noded traingle as 
sh o^^m i n f iq o (3 . 1 ) o 


TfiiO approximations for u, v, p, T are as follows: 
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MV V' 


(3*2) 
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(3.3) 


and 




Tlie vectors V®, T®, and Wy, W®, ¥® contain the 

nodal values of velocities, pressures, temperatures and 

the corresponding v/eights. 

' *' < 

The expressions for the shape functions for velocity 


and temperature are g.lvGn byt 



=- § ^ (2 - 1 ) 

- ( "I “ ”■ ^2 ^ ^ ^ ^ ^1 " 2^ 

N4 == 

I'J^ = '1 ^r, ( 1 “ *“ ^2^ 

1 % - (1 - - ^ 2 ^ 


The pres'inro shnpe functionr-, I'^ui'i: 

T , 9 f • 

3*9(1; 

vrhere I ^ §2 ‘’^''' J-ocnl cnnrd.lnntos and are related 1 

the globalcoordlnnten x and y l)y the fnlJovrlng express.i.onnr . 

^1 === ®1 ■' ^ 2 ^ ■'■ ® 3 ^ I 

*2 “h ^ I 

i 

vrhere | 

= '9 h - - x,y3)/2B« ! 




'1 


,f - 1 


^2 


- Y, ) /2 O5 


®3 = >"3 = - -''3)/20e 

and 

Q^ = nrnn oF thn o.l'^inont = ^ £ >^ -| (yp-Y^)+j^2 (y^^Yl ) ? 

+ ^3 (y^ -IJ ) . 

*2 E.l^mr>ntal A ro n Hxpri-’ss ? ons 

The olomontal area expression a are the contri bntions 
of the eloinent to the left— side inati'ix. Before derivinof 
the eleniontal area ecf’3'’^tions , "the f ollowing vectors d'v, iji 

for achievinrf cornpnc Lnoss in the o<'(ncations» 


Temper <atnre ^^radient vector J 
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where 



Also the following matrices aro defined: 
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Substituting equations [3.2 -- 3.6] into equations [3*10 - 
we get, 

V i = H 

e = i 

and e - R 
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( 3 ^ 15 ) 


on observinq the oquatlon? [2*31 - 2o33 ] from Chapter II it 
is noted that thorn am some non-Unonr terms* These 
laon-linear terms can be made lin'^nr by calculating the 
terms tising the previous iteration valtics. For example 
the non-linear term in equation [2«-31 ] is linearised as 
f ollows : 

a is calculated from the ten’-perature T calculated 
from previous iteration and Vis approximated as V^^ 


v/here 


-n 

V 


similarly for 





NV 
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(3»16) 


The other non-linear terms aro also treated along 
similar linos in the equations [2«32 - 2o33]o 


Substitutinq equations [3.12 - 3«16] into equations 

, e 

[2.31 - 2.33] and integrating over olomontnl area q , 
the following area equations results 
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on obscrvinn tl'^o eqnntlnnr! [2*31 ~ 2*33 ] from Chapter II it 
is noted that thorn aro rjoruo non-Iinnar terms. These 
jjon-linonr tm-'ins can he marie linr^ar hy caJ.culating the 
terms using the provions itornt.ion values. For example 
the non-linn.ar term in orinatien [2.31 j is linearised as 
f ollowG I 

■"a is calculated from the teusperature T calculated 
from previous iteration andi V is nppvoyi mated as 


v;here 


v 

n 

similarly for 


MV V® 
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(3. Id) 


The other non-linear terms arc also treated alone 
similar lines in the equations [2*32 - 2o33]*’ 

fubstltutinq efiuations [3.12 - 3.16] into equations 

e 

[2.31 “ 2.33] and integrating over olemontal area q , 
the following area eqtiatlons result: 
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The elemental matrix j e 
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To convort the above orntatlons from the global 


coordinator to local coordinates x 

and y arc changed 


i n to £ -I 

and ^2 * 
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whore J d'^torniinon t of fhc Jncohi.on inofrlx 
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Tho ohovo oxp.voGolons are then evnlt-iated by 7-point 
Gauss intocfrat:i on scheme for triangle. 
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where l-U ~ nufnb rr of* gnuss points and V/. = W.sights. 

Q 1 


3.3 Global Assciiibly of Hleinentel. Area F.xprossion 

^ £> Y.f X» lip* ily M.j' cf.lobal matrices 

containing tho nodal values of prf’ssnre, velocities, 
temperature and the corresponding vj'-- 1 (Thts . Define 
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The mahricGS [vj f Ey2 hss obtained tiom 

corrcsponci) ng !’;vl mntT"icos ^ » "■|2'*****^3 
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v;here r, S nro tho 'Tfiolml node numbers of the local nodes 
p,n of o] ornont e. 

C is tfio connoctivn ty ■Tor nodal velocity and temperature 
varla’rlos and C is the connectivity for nodal pressure 
variable. 


3 «4 Element Boundary ExpT'ession 

Considoi' a typical boundary olotriont (fig. 3.2). live 
variation of u, v and T along this is also quadratic. Thus 
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!!rare 11 t V , T and V/^^, arc the nodal values at tho 

houncbary local nodes. The boundary shape Functions are ni'cnn 
by: 
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v/hero Sis natural coordinate of the boundary element 

S r- (2 T- ij^)/ 

and Ij^ “ length of tho element b. 

Substituting tfie equations [ 3*30 -- 3*35] into 
equations f 2.34 - 2.36] from Chapter TI, the following 
boundary olomcnt expressions rosuit: 
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To ovaJaialn those Integrals, first we transform 
th®' coordinates from x to I 
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and then tmo the throe polnh 1 -D fanss integration Scheme. 
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/ F(l)d g 
-1 


‘^qb 

B F{{> 

i=1 ^ ^ 


v/hore arc the wo’iqhts, ir; the number of gauss points 

and are the qnu.ss points. 


3.5 Evaluation of rj*"*, an d o n the Houndaries 

Re-Forrlno to f.lqure (2.3), it is observed that 
convective and cnniluetive boundary conditions occur on 
sides 2,3 and 4. On f’ther sides, it is either unspecified 
or ;!:Gro as on side 6« 

As on sides 2 and 4, the expression for non-dimensj on- 
alised convection is given as: 


q 


n 


h h 

. o 


(T^ - 1) 




From equation (3.32) Tj^ can bo expressed as: 


— T -b 

Tb = I 


Substitution of the above equation in equation (3 .41) 


results in: 


£ 


4d.- 


(3.42‘i 


•o *'0 b 

on side 3, the non—dimensionalised heat xlux equation is 


given by: 


'n 


-q A A + s. ^ 


where, 
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C- 

4 .* 



2 r K 


The evaluation O'f^ the terni ie dep'^ndent on whether 

the region conaidorod •* o AB or (''--'e fig* (3*3)) i#e., 
in hackv.’nrd ollp region o.r -Torv/ard slip region 

backv;nrd ’H )■) region, the expression for t^, is 

-- t iirv^ -*■ t 
'•s X y 

^ u din ‘3 ~ V CosO 


and 

rr - ■'iin‘^ (u :un<^ - v Cos<^) Cos9 in' SinP - v ■•> v 

-)■ C,^ (f - 1 ) 

jti 

(;a„ 4 Pn;;. 


and similarly in forward slip region 


q = C. dinf^ (u SirP - v Co; 
n ' 


^D> t " CPlosQ (u Sin© - v CorP) t 
X • 


C2(T - 1) 


(3 •43! ) 


vSubs h Ituting equations (3*43) and (3.30-3*32) into 
equation (3.40) vxe get for backuiird slip regions 


T ^ . b , r,.. ^ ^ r n . b 
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Th <7 expro'v; I on ;' for jj contain the unbnnwn nodal values of 
b h .J'* - _ 

li » .Y. , T and t^. Thd'so to-ni!.s aro t.airon to the loffc- 

hnn.l siJn ,-,rtnr nvn Inn Ling nnd from previous iteroLi on. 

nv.ilmtion O'- nnd r,’’ is oonrin.O'1 now only to .oido p. 
As on othd'r slilos it is oJthor un.spoc.i C led or zero. On 
side 3, \io Iiavc a rnlatlon b^'twoon and t^. This rGlat>'''n 

i. s dnpondonb on v/hothor the region consldnrod is backv;ard s.Mn 
or forv/ard slip rogion. 


For roqlon A3, tiio relation hotv/een t and t i^* 


Y 


IS a: 
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where t 
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normal fn -co. 


(t^ FosO + t sinQ) 
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— ^ ...p ^ ^ .~™x. 


ty Slna 1 n Co^Q 
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or t^ ” t Cob*^ 

X V 


O y 

^3 {.I w r 
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( » 1 5b ) 


Umilarly for rogion nc, wo hove 
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Coaf^ •! T) .'jinO 
3i.n^ " tlCos*^ 
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3.6 G1 oi-i a. T A a ^ ^ ^.rHir irin x'y Expressions 


I/'t. , 




J Ix 
''q 0 


y 


(^.J ^ + ig ) 


(3,47) 


v/horo tho mabrices W and ^ have boon dorived earlier. The 
matrix Q. ,is the aaaoml'ly of the knov/n terms in the 
expreo-sion for £% on the sides 2,3 and 4 \^hile the coeffi- 
ejonts of u,v^ and conatittibo to the matrix A-| * fhesG 
matrices are oi.italned by expan'll ng the appropriate terms 
and adding over the reIo'«/ant boundary elements. 
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3.7 Globnl I-gun I: i on 

ComhiriJnqf the orjuntion?? [3.3/’] end [3.47] we get 

K ^ WT + B ] 

T 

3 1 nee W / Oj vjg heve 

K ^ ~ A’] ^ 1], 

or A A; “ I? 

wi'i ere A ~ K ~ (3- .4 3) 

3 • ^ Botm ei ary Condi t'i enr; 

3 .B.1 Spe cial noundar y Co ntb*. h '.on e 

In thic section the application of the mechanical 
boundary conditions on sido 3 is dJccucaed. These condition'^ 
Aeinq the equations [3.4b<a, 3.4Bb] atui [ 3.46a, 3.46b] and 

4 - 

i.e., u + V tan 0 = 0 

These boundary conditions are applied by making some 
rov/ operations on the middle set of equations (equations 
corrospondinq to Y.} • Tnltially in the equation corresponding 
to a node on side 3, both the right hand sides are unknown. 
Then ttsinq ociuatlisns [3*4b3, [3«4b] and [3o49] we make the 
right hano’ s1 de corresponding to u-oguatlon known. Then 
u-equatlon is r-nlacod by the ex prosed on [ 3 .49 } . 



3.8.2 Er^.q.nntx.il nonnc j -iry Con. l UI on^ 

Tho Gosenfcial honn;iary c^n ilt.i otin are applied by 
finding out the rov; correaponding to that variable in the 
global matrix arul mating it all excf’pt for the diagonal, 

which ia made into un i.ty. The corresponding value of the 
vanablp is insortod in the riglil hand side. 

Rrr applying u-volocity condition at node R 

Wo make A (MOp-i , J ) --- o for J = 1,N and J /- IPp-i :> 'dV- 

“ 1 for J N0C-(2^R-1 

dlmila'sly for applying V-vi^loc;ity condition at node R 

v;g make A (MOr:-i2>f-R, j) - 0 for J , M and J / nOf>ip’<g 

1 for J ^ M0C+2*R 

For applying t'^mperature boundary condition at node R 

we make A (M0C-i2’‘'N0T+R, J) = 0 for J=1 , N and J /MOg-i-p^P'Ori 

=- 1 for J = NOG + 2 ^ hot -> r' , 


3 Progrnmise Implement at ion 

Tv/o programiios were developed in the process of 
solving the problem. 

In the first program a non-frontal technique was 
adopted. The assembled global matrix v/as fully stored 
because of the lack of symmetry and bandedness* After 
applying all the hotmdary conditions, NAG routine F04ARF 
v;as used for solving equation. This routines calculates 
approximate solution by using Grout's Factorisation moi.hod. 



The draubackn v.d i.fi i:h'j proqram vms the huge demand it plac es 
on memory and CPU t1n;a. The typical time it took for 
10 iterati'ans v;as around 2 hours, 

Th o cecond t-'Tnnicun was doc'^i nped liecause of the draw- 
backs in tlie eaiU.ior program andi [-,h<'’ ' n ak”-.! llty to increase 
the elements further- This program incorporated Frontal 
technique [11] for solving tho ogoali orir . This program 
was tmtcfi Tocher tliian 1'ha oarUlot- program and occupied loscer 
core tn'^mevy. The only limltatien being the disk space. The 
main anci c''uci.al iJravjdack \ji. l:h this pirogram v/as that it 
does nor. do complehe pivoting. Tills load to erroneous calcu- 
lations oT nodal pressure voluos because the diagonal term ii 
the continuity equations wore very small owing to the penalty 
f ormui a i; I on • 

Flna.lly it was decided to stick on to the earlier 
program using NAG routine. All the results presented here 
wore ca.iculated us.lng this 


progT'prn 






4,5,6 are the mid points of the corresponding sides 
Flg.3.1 Typical element e and the local node numbers. 



1 


Fig, 3.2 Typical boundary element b and the boundary 
iocol 'node numbers . 




AB- forward slip region 
BC- backward slip region 
horizontal datum 

Fig. 3.3 Backward and forward slip regions on the 
roll -work Interface. 
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TV 

TnMPLiPATTRp ; .! 5A'~'’ ^ iRF;: rpiT T r-i n ni. l_t ? n 

fn th i chnphor the c^xporj in^ntn I setup along w:! th 
the i!ir>!-l(n<! ,-i(inpt^u jn moasnrlnp the t.e-.'.perature is 
d i see s sed. 

T'ho roll'ni'r nsod uas a two-!iigh roll.i ng 

in eeii i vi [:!i li'vr. K^rv'd .»ne’ lepp'^'l ’’’oils. The meehnio 

has e, rv^'''en''('-'iv,nn t to vary T.hn retiuetlon ratiOe 

The nsetho'i oT o[!ihe(ido(j theriiiocon pie was used XU 
the ir,easu''''!rent . A Ch I'oinol-Alou', coinh.i nation of thoxs -o- 
cotjp.io vras e;r,!ied.]eci on the v/oxir in'rr,o'’'i'il as shown in 
’“ig. 4.1. nofore the the.rinnconplo v;as embedded, it v/n^ 
coated v.d til a etminel paint « This frn'n:ed 'as a insulatinrY 
cover on the vrlre and prevented the the.rinocouple wire T.ro<r, 
corning in direct contact with the v/ork material. ■ The 
leads xrore then taken out and connected to the recorder. 

VJhen the slieet v/ns rolled, the thermocouple embedded in 
the v/ork sliovred a temperature rise on the recorder. A 
sketch of the setup 1s sliovm in Fig. 4.2* 

■'lie above experiments were conducted on three pl--^ccs 
of aluminium strip for three different reduction ratioso 
The sp'^ed of the roll v^ras kept a near constant in all the 
experiments. Th^ exper unental iy, measu.rod values and 
theoretically calculated values of temperature are pronented 
in next chapter. 




1 - Strip metal 

2 - Thermocouple wires 

3 - Thermocouple bead 

Fig. 4.1 Arrangement of thermocouples on the work- 
piece . 



1 - Strip metal 

2 - Thermocouples 

3- Thermocouple bead 

4- Roll 


5 - Recorder 
6 -Pen 
7 - Graph 


Fig. 4.2 Schematic diagram of experimental setup. 




CT-IAPTER IV 


RE‘7rLTS AND UI -SCySSION 


This chapter presents results for the variation of 
velocity and temperature along the interface. The 
temperature obtained by experiments are also shovm. The 
tv/o materials used for theoretical analysis were Aluminium 
an;i ntee.l. Experiments wore conducted only for Aluminium 
sfieet ov/inq to the limitations of the rolling machine. 

?\~‘tenial Properties 

’following m'o the thermal and mechanical properties 
v;hich have been used in the analysis. 

Muminiinn 

''iolci strength = 88.0 (MPa) (assumed to be constant) 

y 3 

Density at room temperature “ 2707 (kg/m ) 

Tno'-ficient of linear =a = (-91 .0+0 .38 T)x 1 0”^'(por°C) 
pxinnsion 

'specific heat = c = 674-I-0.6T (J kg per C) 

Thermal Conductivity = K =239 (W/m-°c) (assumed) 

m to be constant) 

3teel 

Yield strength = d = 217-0.15T (MPa) 

at room temperature =^q “ 7870 kg/m ^ 

Coefficient of linear expansion =a = (-48+0 .2T)pGr C ; 

TpociTic boat = c = 296.6 -t 0.48T J/kg- C ; 

Therm, .il conductivity = " 


76.8 - 0.01T 
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Iho roll matorlnl ur.eci was inedlum carbon steel* Its 
thr-'n-.a] onn-inctivlty \/as taken as 30 Joules /m/°c/sec. 

Those vaJiios were taken from TPRC Journals. 

The program v;as executed using a 56-element mesh 
sh'^'i'/n in ’-ig. 4.1 o Tlie results obtained experimentally 
for aluminium for reduction ratios 15, 25 % reduction 

compare reasonably \;ell with the values obtained from the 
th o or r-t-. 1 c a 1 mo do 1 . 

■"ho theoretical analysis was carried out for five 
■ b' r sorent reductions of Aluminium and tv/o different reductions 
o steel at various angular velocities of the roll. The 
graphs for temperature and velocity distribution along the 


interface 

are presented. The 

details 


Aluminium 




Case 1 : 

Initial thickness 

= 3.25 mm 



Pinal thickness 

== 2o75 mm 



Reduction ratio 

^ 15.3.R'(J 



VJidth 

= 50 mm 



Angular velocity of the roll = 2.0944 

rad/sec . 


Poll radius 

= 65 mm. 


The 

graphs for this case arc shovm in Figo 

4.2a, 4.2b 

and the experimentally record 

value is depicted 

in Fig. 4. 2c 

Case 2: 

Initial thickness 

= 3 o25 mm 



Final thickness 

= 2.45 mm 



Reduction ratio 




Width of the sheet 

= 38 mm 




} 


Anqn.inr volocity of the mil = 2.4 rad/sec. | 

foil radii 10 = nrr. I 

! 

I 

The corrf^r.ixindinq graphs for this case and f^igs«4.3a, 

4.3I> an'l th<^ orpori inontally recorded va-1uo in Fig* 4.3c. j 

g.’ so n : Int hinl l-.hi cknoss - 6*274 mm 

'■'inal thickness - 4. *^00 nim | 

n-'d.uction ratio :r. 21 

'.Ji.'th of tlio niieet = 38 mm 

Anon.] pr velocities of tho roll being 1.92308,2.40334 and 
3.3461 Rndksec. 

Roll radius = 79.375 mm. 

Tho corrosponding graphs for this case are Figs. 4 .4 - 4.6. 

gpso 44 Initial thickness = 6.274 ram | 

Final thickness = 4.445 ram. | 

Rf'dnctlon ratio - 29.40-^ , 

V/idhh of the sheet = 38 mm 

Angular velocities of the roll being 1 .92308, 4.40384 
and 3.3461 Rad/sec. 

Roll r<a..iius = 79.375 rnra. 

The corresponciing graphs being Figs. 4.7 -- 4.9. 


3 a s o ft j 


Tnl-i-ial thickness 
s*inal thickness 
Reduction ratio 


2.032 mm 
1 .346 mm 
33 . 75 ?^ 


Rhdth of the sheet = 33 mm 

Annular volocitios of the roll being 1.92303,4.40384 
'and 3.3461 Ra'i/sec. 


loll radius 


7<^.n75 mm. 



The figures fo- this being 4.10 - 4o^2. In all 
th*^ nbovp cases o" aluminium specimens, the frictional 
censtapt v;as assuraeu to bo 0.2. 

steel 

C ase 1 t TnitTol thickness = 1.0 mm 

-inal thickness = 0.84 mm 

Cediictlon ratio - ITT-J 

Vfidth of the sheet = 50 mm 

Annular velocities of the roll being 2.015748, 4.031476, 
5.03937 Rad/sec o 

Roll radius = 65 mm. 

The figures tc*' this case are shov/n in Figs. 4.13 - 4.15. 

Case 2 t Initial thickness = 1 .0 mm 

Final thickness = 0.76 mim 

Reduction ratio = 24"^ 

Width of the sheet = 50 ram 

Angular velocities of the roll being 2.015748, 4.031496, 
5.03937 Rad/sec. 

Roll radius = 65 mm. 

The figiires for this case are 4. 16 - 4.18. 

The frictional constant for the above three cases 
being 0.08* ; 



D'i .‘'0 i nn 


An ox nr'<5 notion of the t erripernturos measured by experiment 
and px'oiHctod by tfioory clearly shovr that there is reasonable 
aaroement* The tlieorotically predicted temperature profile 
ns sho'.ni in nil tlie cases shov/s an initial peak rand a gradual 
decrease along the interface towards t’ne exit. The experi- 
nw'-d-ell y mea-snx'ed teinperaturo pen Pile also shov/s a propensity 
Onr ciecr nefri.ng (along til'’ interface) towards the exit* Tt is 
ni so olxsG^ved that the decrease in the experimental temperature 
values along t!io interface is more gradual than the theoretical 
temperature values • 

The velocity profiles as observed from the graphs on 
tile interface sliow tliat there is a initial decrease and a 
gratiual rise later tovjaj'ds the exit. This might be attributed 
to the coarse mesh adopted in the onalv'ais* The location 
of neutral point v/hero the value of tho non-dimensional tangential 
voloci ty V is egunl ’ho uni ty is in good ngreemont v;ith the 
neutral point location calculated fi’om Avitzur*s and 5im*s 
formula . 

Some of the special features noticed in the results 
are the decrease in the interface temperature towards exit 
contrary to what most of the researchers had predicted. 

Altan predicted almost a constant temperature distribution 
a. long the interface. Jesweith and Rice measured the 



fcoinporatnro rise pxpprimontnlly anci found the temperature 
f;o ri ,r;o r.llphhly In.i t tally and then even out along the 
interface tnv;ards exit* T!vl a mioht be dne to the sticking 
f^'icfclon concfcrainh put by Alton in his analyses. The 
cainpartson utth Jos'./oith and Rice’s method leads to 
canclucion throt the thorrnocouple used in the present 
work t'of(uirod longer time to respond. The other major factor 
bo'irig hho method of m onsurornent . Tn the present vrork thermo- 
couples \jorc omlecdded in the v7ork. In Joswoith aid Rice 
work the thermocouple was embedded on the roll. Thus in 
our case the thermocouple bead always v:as lagging in velocity 
vn'.th roll and Crome in continuous contact vath the cooler 
parts of tlie roll unlike in Jesweith and Rice’s \;ork. 

The nossiblo explanation for the drop in temperature 

rise along the interface towards exit as seen in the graphs, 

imy he obtained from observing the velocity profiles 

along the interface. It is seen that the velocity gradient, 

d^/c; 

i.e., strain rate (~^) is decreasing towards the exit. This 
results in reduced plastic vrork and hence a drop in temperatur 
ri se • 

It is also observed fram the graphs that the temper- 
ature rise increases vath increase in reduction ratio and 
angular velocity of the roll. The temperature profiles are 
smoother for lovier velocities of the roll. The change in 
the velocity profiles along the interface with increase in 
angular velocities is hardly perceptible. 



"or -utijre Uorh 

In 'I'ho preoont v/nr!c attempt.'^ inxe made to predict the 
roll-interface prossu-fe di '^tribu tion, roll-interface frictional 
"n?'ce 1 !i;'t\'ibntlon and roll torque but wore unsuccessful. 

This vras nvanq lo the erroneous calculation of pressures, 
bo a ■further extension to this v;ork vmuld be to calculate 
th'^'Si"! naratnoters with both back anfi front pulls and arrive at 
some optimisation. Another useful extension vrould be to take 
Into account tho offeebs of strain rate) and total strain 
■ n c •"’J.culati on of media ideal properties of the material. FuJther 
scope also lies in not assuming the 2-D constraints and 
usi.ng b— D PHH technique predict the spread of tho sheet 
metal and the sheet profile actor rolling- 



























Fiq.4.3b. Variation in Temperature along the Interface 

Reduction ratio = 2^ PaH/cor* 

Angular velocity of the roll = 2.4 Rad/sec 



, Cticrn, - 



t X PEKIM tt^T ALLY Recorded G)R 


4.5c 


^FnnrTtnM - 







Fig. 4.ifb. 


ike OisUnce in H 

Variation of Temperature along the Interface 


Reduction ratio = 

Angular velocity of the roll 


1 .92308 Rad/s f 

























Fig. 4.?b. Variattei^^5^®'^lilf)Srature along the Interface i 

Reduction ratio = •ifGX | 

Angular velocity of the roll = 1.92308 Rad/se 
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Fig^ 4.9a. Variation of V along the Interfaces 
Rediction ratio = 23 •^% 

Angular velocity of the roll = 3.8461 Rad/sec 


jv'4 ^ 




fire Distance In H 
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Fig* 4.9a* Variation of V along the Inteirfacei 
Rediction ratio = 

Angular velocity of the nil = 3.8461 Rad/sec. 






Fig- 4 3b. 


Variation of Temperature along 

Reduction ratio = 23 -4^ 

Angular velocity of the roll = 


the Interface 
3-8461 Rad/sec 





fire Distance In M 


Fig.4.10a. 


Variation of along the Interface ^ 

Reduction ratio = 33.75% n j/ ' 

Anqular velocity of the roll = 1 .92308 Rad/sec 
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ftrn Dlstidnce In 11 


1.3X10"^ 


Flq* 4.10b* 


V.irn at.i.on <a-^ tamFerabJr-'^ along the interface 
P.o^juction ratio = 33*l-'fo 

Angular velocity or the roll = 1 <>92308 Rad/ sec 
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I'-iflation 111 1 ’:? aluiiq the interface 
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A'txi.'irvi -ni M.-rr) 
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1.3X19’^ 
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riq.4.11n* Variation of V^ along the interface 


Reduction ratio = 33.75fo r> ,/ 

Angular velocity of the roll = 2.40384 Rad/sec. 





, ^^'irlation cf Tcspmi^ure aion*! the Interface 

r— I Y j .j. I j 


tVLu.'vru.'ru. ultt) 



Dlr»^!iin:B In II 


F5,g. 4*1 1b. Vari-ition of temperature along the interface 
Ro'ltictl-on ratio = 30 .75/o 

An pilar velocity of the roll = 2.40384 Rad/sec 




fin: !n H 

FH*4.12a. Variation of V alona the Interface 
Roiitictlon ratio = 33.75''o 

Anqular velocity of the roll = 3.8461 Rgd/seco 



1 r> Km 



Arc Dis^swe in M 


’^Iq. 4.12b« Variation of temperature along the 
interface 

Rpdtiction ratio = 33o7o,o 

Angular velocity of the roll = 3o846l Rad/sec • 




ftrc OlSlrSnne In H 


I'lcf. 4.13a. Vnr1.ntlon of Vg along the interface I 

Reduction ratio = 16 % ; 

Angular velocity of the roll ^ 2.015748 Rad/sec. 
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3.2X18^ 
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Uariatlpn of Tefaperat ure a jony }\ l e In terface 
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Hn; Dlpt3f»c0 la H 

Fiq* 4,13b* Variation of temperature along the interface 
RoiJuctlon ratio = iQ’o 

Angular velocity of the roll = 2*015748 Rad/sec* 
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fhr In 11 


.|.14a. Variation of along the interface 

AnSaar^volooity of the roll = 4.013496 Rad/sec 






ftrc DlDfcf»nre hi H ,, . ^ ^ 

4o14bo Variation of tomperature along the interface 
Reduction ratio =16'- . n j/ 

Angular velocity of the roll = 4.0134^0 Rad/ 



^ mp ^ ^ ^ 
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Ficr. 4.15b» V-^riation of temperature along the interface 
Rr^duction ratio = I6”!l 

Angular velocity of the roll = 5.03937 Rad/sec. 




Afp In II 

^ig. 4.15a. Variation of V along the interface 
Reduction ratio = 16?o 

Aniular velocity of the roll = 5.03937 Rad/sec. 






An; In H 

Fig. 4.16b. Variation of Temperature along the interface 
Reciuction rati o = 24°?; 

Angular velocity of tho roll = '2«015748 Rad/sec 



Variation of Us alonq the Interface 



Arc 0i.ihftnrp In H 

Variation of Vg alonq the interface 
Reduction ratio = 24% 

Anqular velocity of the roll = 2«015748 Rad/sec 














RnFFlRGMCEG 


1# Jnop C. fbqt'v.j.ial and F-'ran',:; FF. Valdpaus 

"On tFio liisplo^’i'^ntation of finito atarain plasticity 
onna {; Ions in a numerical modol", ilumerical Analysis 
of I'n"!r,lrrj PmcesseG, 1934, p® 351 »• 

2. G,D. Lalioti, .'i.’F. Siiah and T» Alton 
"Cuinnutor aided analysis of the deformations and 
teir.poraturos in strip rolling" , Journal of EnginGcr ing 
for Industry, Trans, of ASME, Vol. 100, May 1978, 

n. 159. 

3. 7.i enkiev/iv/, Onate and fleinrich 

"A gr'noral f orn-.nlati on ^nr coupled tJiermal flow of 
mot'ils usinn f inite e lemon hs" , Ini;, J. Mum. Metli« 

Eng . , Vo 1 . 17, 1 9p i ^ j., , 1 4 07 , 

4 . n,9. Co’f'nf^iold and Fi«!!» .T^Jmson 

" Tiionro tical prediction of pilnstic flow in hot roiling 
including the effects of various temperature distri- 
Isutions", J. Iron A Steel Inst., Vol. 211, 1973, p. 567 

5* V/o Johnson and V. Kudo 

"The use of upper-F;iound solutions for the detennin tst ion, 
•of temperatt-ire distributions in fast hot rolling and 
axi symme trie extrusion process", Jl. of Mech. Sci«, 

Vol. 1 , p. 175. 

6. J. Jesuiet and VF.B. Rice 

" M'^asurement ot strip temperature in roll gap during 
cold rolling", Annals, of the CIRP, Vol. 24, 1975, 
p. 153. 

7. A.M. Karagiozls and J.G. I.ruord 

"Temperature distribution in a slab during hot roi.lina", 
J. of Eng. Materials and Tech., Vol. 110, Jan. 
p. 17. 

8. Avitzor 

"Metal For '.sing Process". 

9. Chaturvedi 

"Rolling of Sheet Metals". 


to. 

J. 

FA Mo 1m 

■^n 









in 7 

'^ai: Tro 

I'v*; r 

AT” , 

Me' 

' 'C 

•u,; dill 

V‘ 

a, tin n 

, 1981 . 

1 1 . 

n 

* 

Tnyi Of' 

n^i 

'j I-’' # 

'Mr 







t5 l\ 

1 \ 

ntim''''ri 

■ : a 1 


" !;i on 

0 f tl '! '• 

!F 

,'ivi 0 

S ok e s 



''O 1 ; i ons 

! 1 "4 

1 rvi 

f:he 

f 

;ltvi tf' e 

loi 

i'’nr\ !; t 

ociuuiquc 



J to'ss 

1 

Cl 

Ui d 

f 

Vol. 1 

t 


p » 73 • 



D.Mo Irons 

"A frontal solntion program fo3: finite element analysis 
Int. J. riumo Meth, in EIngg., 2, 1970. 

J.Mo Reddy 

"An introdnction to the finite element method", 

.'.'loGravr I 'ill Book Company, 1934. 

C,C» /'.ienkiGwic z 

"The finite ele>r.cnt method", Tata McGraw .Hill edition, 

1901 . 

LI and S. Knbayashi, "Rigid plastic finite element 
analysis of plane stron rolling",. J. of Engg® for 
Tndnslry, Vnl. in, 1982, p© 55© 



